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We present experimental, numerical and analytical results for the study of near-field interaction 
of twisted split-ring resonators, the basic elements of the so-called stereometamaterials. In contrast 
to previous results, we observe a crossing point in the dispersion curves where the symmetric and 
antisymmetric modes become degenerate. We introduce a model to describe the interplay between 
magnetic and electric near-field interactions and demonstrate how this model describes the crossing 
of the dispersion curves, initially considering lossless identical resonators. Finally, we apply the 
theory of Morse critical points to demonstrate the competition between losses and fabrication errors 
in determining whether or not symmetric and antisymmetric modes cross. 



I. INTRODUCTION 

Metamaterials created as an array of sub-wavelength, 
resonant elements can exhibit interesting electromagnetic 
properties, such as a negative refractive indesP'. An im- 
portant building block of metamaterials is the Split Ring 
Resonator (SRR)^, which has a highly dispersive mag- 
netic polarizability which is negative over some frequency 
band. However the response of a metamaterial is not sim- 
ply given by the response of an individual resonator, but 
depends also on the arrangement of resonators within 
the systerrP. Unlike the atoms in natural materials, the 
near-field patterns of metamaterial elements are quite 
complex, giving rise to strong interactions between them. 
This means that the arrangement of resonators in a lat- 
tice or super-lattice plays an important role in determin- 
ing the response of the metamateriaP^iSl. A review of 
many examples of interaction between such resonant el- 
ements can be found in Ref. [TTJ Characterizing the un- 
derlying interaction mechanisms between elements is es- 
sential to understanding the overall resonant properties 
and effective parameters of the material. By controlling 
the relative arrangement of elements, it is possible to al- 
ter the response of the material substantially even if the 
constituents are fixed. 

Recently, H. Liu et aP^ analyzed numerically an array 
consisting of pairs of split ring resonators on the same 
axis, with the second ring rotated by 90°. The hybridized 
modes of this system show strong polarization rotation, 
and have suppressed radiation losses^"^. Subsequently, 
N. Liu et aPSi investigated a similar system, but consid- 
ered an arbitrary angle between the two rings. Using nu- 
merical analysis, they extracted the eigenmodes of a pair 
of SRRs operating in the near infrared. They concluded 
that as twist angle increases, the resonances converge, 
undergo an avoided crossing, then diverge again. The 
reason for this avoided crossing is not clear, since our 
previous study of a different systems of coupled SRRs 
showed that the hybridized modes can cross- . In Ref. [T3] 
the dispersion curve with an avoided crossing was fitted 
by a multipole interaction model, however no physical 
justification for this fitting was given. In particular, it 
was assumed that the magnetic interaction constant is 



invariant with twist angle, however this is inconsistent 
with the strong variation of current density around the 
circumference of the ring. This model was subsequently 
extended to include the polarization-rotation of the scat- 
tered radiatiorP^. 

The purpose of this article is twofold. First, we give 
a rigorous analysis of the interaction between a pair of 
co-axial SRRs with one ring rotated, and determine the 
conditions under which the dispersion curves will un- 
dergo a crossing or avoided crossing. We introduce a 
physically-based model for describing the interaction be- 
tween the rings, and show how this model explains the ex- 
perimentally and numerically observed dispersion behav- 
ior. In Section [llj we present our experimental and nu- 
merical analysis of the dispersion curves. In Section [ill A| 
we show how the dispersion curves can cross, based on 
an idealized model of a pair of SRRs. Finally, in Sec- 
tion IIIIBI we consider the influence of losses and non- 
identical rings on the dispersion curves. Using the theory 
of Morse critical points, we will show that competition 
between losses and fabrication imperfection will deter- 
mine whether or not the modes cross. 



II. EXPERIMENTAL RESULTS 

We consider a pair of SRRs with varying twist angle 
9 between them shown schematically in Fig.jlja). First, 
we perform microwave experiments with a series of rings, 
with one ring held fixed, and a separate sample created 
for each rotated ring. The rings used have an inner ra- 
dius of 3.5mm, an outer radius of 4mm, and a gap of 
1mm. They are copper, printed onto 1.6mm thick FR4 
circuit board, and the rings 3.6mm apart, with the di- 
electric boards located between the rings. The incom- 
ing microwaves are polarized so that the electric field is 
across the gap of the fixed ring, as shown in Fig.jlja). To 
match our experimental data, we preform numerical cal- 
culations using CST Microwave Studio with the boards 
having a dielectric constant of 4.6. 

Experimental results were measured using a Rohde and 
Schwarz ZVB network analyzer in a WR-229 rectangular 
waveguide. We measure the excitation of the rings from 
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the absorption of the system, given by 1 — 1 P — | Si i ^ , 
where S21 is the transmission coefficient, and Sn is the 
reflection coefficient. A comparison of the experimental 
and numerical absorption curves for 9 — 90° is shown in 
Fig.jljb). To characterize the dispersion behavior, exper- 
iments were performed with varied from 0° to 180° in 
10° increments, while numerical results were calculated 
in 5° increments. For each angle a Fano functiorii^ was 
fitted to the experimentally and numerically obtained ab- 
sorption curves, and the resulting resonant frequencies 
and absorption coefficients are shown in Figs, and 
(b) respectively. In all cases the data were well-described 
by a pair of Fano resonances with a correction for back- 
ground absorption, and the fitting to the numerical re- 
sults has a residual error below 5 x 10"'^. 

For 9 = 0°, there are two resonances ojs and ujas, 
and by inspection of the simulated currents in the rings 
we verify that these correspond to the expected sym- 
metric and antisymmetric modes. As 9 increases, ujas 
increases and ljs decreases, reaching their maximum and 
minimum values respectively at 9 = 180°. As can be seen 
in Fig. [2jb) and Fig. |3][a) , the two resonant peaks have 
different widths, which are primarily due to differing radi- 
ation losses. The symmetric mode has relatively stronger 
radiation losses for low angles, since each ring approxi- 
mates an electric dipole, and a pair of parallel dipoles ra- 
diate strongly. As the angle approaches 180°, the dipoles 
become oppositely directed, thus we have the low radia- 
tion efficiency of an electric quadrupole/magnetic dipole 
like distributiorfl^. The antisymmetric mode has the 
charges on one ring of the opposite sign to the other, 
therefore it changes from an electric quadrupole to an 
electric dipole type distribution with increasing angle, 
and the radiation efficiency increases. 

Inspection of Fig. [2]ja) shows that the resonances ap- 
pear to cross at 6* w 34°. We can be confident that 
there is not a very narrow avoided crossing, by consider- 
ing the corresponding numerical loss coefficients shown in 
Fig. [2|^b), which are clearly distinct from each other near 
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the crossing. The numerical analysis was repeated in 0.5° 
steps near this crossing to verify that no features have 
been missed, thus there is a clear indication of a cross- 
ing of modes. In Fig. [sjjb) we have plotted the numerical 
absorption spectrum in the vicinity of the crossing. It has 
the appearance of only a single peak, however the fitting 
procedure still identifies two separate resonances, which 
are also plotted for comparison purposes. In general the 
experimental data show good agreement with the numer- 
ical data, but with some uncertainties due to fabrication 
tolerances. We will show in Section [ITTB] how experimen- 
tal uncertainties contribute to the crossing behavior, but 
first we will investigate the interaction mechanisms in an 
ideal system. 



III. THEORY OF CROSSING 

A. Identical and Lossless Rings 

The tuning of the system by rotation can be explained 
by looking at the interaction between the rings. As the 
rings are twisted, the magnetic and electric near-fields 
between the two rings change, changing the coupling be- 
tween them. We first approach this problem using the La- 
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FIG. 1. (a) A schematic showing the rings rotated with re- 
spect to each other through angle 8, and the polarization of 
the incoming waves, (b) A comparison of the experimental 
(solid) and numerical (dashed) absorption for angle 6 = 90°. 



FIG. 2. (a) A comparison of the experimental (markers) and 
numerical (lines) resonant frequencies, (b) Corresponding 
absorption coefficients, calculated from the resonance line- 
widths. 



3 




3.4 3.5 3.6 3.7 3.8 3.9 
f(GHz) 



FIG. 3. (a) Numerical (solid) and experimental (dashed) ab- 
sorption curves for = 180°. (b) Numerical absorption curve 
(solid) for 6 = 34° at the crossing, showing the two fitted 
resonances (dashed). 



grangian for a pair of identical and lossless resonator d^ l ^^ l 

£ = A{Ql+Ql + 2aQiQ2)-B{Ql + Ql+2pQiQ2), (1) 

where a and /? are the dimensionless magnetic and elec- 
tric interaction constants, and Qi,2(i) are the time- 
dependent amplitudes of the modes' charge distributions. 
By substituting Eq. ([T]) into the Euler-Lagrange equation 
the dynamic equations are found to be 



(2) 



Solving the characteristic equation for this system gives 
two resonances: symmetric (Qi = Q2), and antisymmet- 
ric (Ql = -(52): 



1 + /3 



WAS — ^0 



(3) 



For a pair of rings in a homogeneous dielectric hack- 
ground^ the electric and magnetic interaction constants 
can be determined from the interaction energy between 
resonators, using the method described in Ref. [B] The 
resulting interaction constants are shown by the markers 
in Fig. |4|^a), along with the following functions which fit 
the data very well: 



13 = Pi cos{0) 



ao + a\ cos(f?) 



(4) 



with /?! = 0.085, ao = 0.098 and ai = 0.05. These 
constants are dictated by the charge separation across 



the gap of the ring, the current circulating around the 
ring, and the inhomogeneity of the current distribution 
around the ring, respectively. 

For rings aligned on the same axis, we expect that the 
magnetic interaction a should always be positive, as the 
intersecting magnetic field from one loop should always 
be normal to the other loop. In addition the electric in- 
teraction /3 should be positive at 6* = 0° as the charge 
distribution has the nature of parallel dipoles. All ar- 
rangement of rings on the same axis which we considered 
behaved in this manner. 
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FIG. 4. (a) Magnetic (a) and electric (/3) interaction con- 
stants calculated for a pair of rings in free space. Dots show 
exact calculations, lines give the fitted function, and (b) shows 
the corresponding resonant frequencies, (c) Interaction when 
electric coupling dominates at S = 0°, and (d) corresponding 
resonant frequencies, (e) Interaction constants which become 
equal, and (f ) corresponding resonant frequencies which cross. 

In Fig. we plot the corresponding frequencies of 
the symmetric and antisymmetric modes, normalized to 
ujQ. As our approach models the response of the res- 
onators in a homogeneous dielectric background, the re- 
sults are significantly different from those observed exper- 
imentally, where the dielectric is inhomogeneous and the 
effect of waveguide boundaries is also significant. In par- 
ticular, for this system of perfectly conducting rings, the 
crossing of resonances cannot be reproduced for rings in 
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a homogeneous background. Therefore, we consider the 
possible regimes of interaction which may occur, under 
the assumption that the interaction constants will be of 
the form described in Eq. Q. 

The case considered in Figs. [4ja,b) corresponds to the 
magnetic interaction always being larger than the elec- 
tric interaction. This results in increasing splitting of los 
and ujas with increasing twist angle, however in prin- 
ciple there is no reason why the splitting cannot de- 
crease. We show such a case in Figs. |4]jc,d), where we 
have set /3i — 0.02 < ai, such that the inhomogeneity 
in the current has a stronger influence than the dipole- 
like charge distribution. Despite the different behavior 
of the frequency splitting curves, there is little qualita- 
tive difference between the interaction constants shown 
in Figs.jIJa) and (c). 

The only other case allowed in our model of inter- 
action under the afore-mentioned physical constraints 
on a and f3 is that a > (3 ior 6 = 0°. An exam- 
ple of this is given in Fig. |4|^e), where we have set 
ao = 0.108, cki = 0.05, /?! = 0.18, such that /3 = a at 
6 « 34°. The corresponding resonant frequencies nor- 
malized to ujo are plotted in Fig. |4]jf). The parameters 
have been chosen to closely match the dispersion shown 
in Fig.[2ja). From this match between the model and the 
experimental data, we conclude that the inhomogeneous 
dielectric enhances the electric interaction between the 
rings but has almost no influence on the magnetic inter- 
action, as expected. 

Equations ^ show that the tuning curves arise from 
competition between electric and magnetic interaction 
constants. If the magnetic and electric interaction con- 
stants have the same sign, then they will counteract each 
other. Thus the frequency splitting can be weak even 
if the near-field interaction is strong. In particular, the 
model shows that there is no splitting when a — (3, in 
agreement with the numerical results presented in Fig. [2] 



B. Non-identical and Lossy Rings 

Up to this point, our analysis has been based on the 
assumption that the two rings are absolutely identical, 
which means that the hybridized modes will be symmet- 
ric and antisymmetric, and degeneracy can occur when 
the effective coupling is canceled. If a pair of lossless 
rings are not absolutely identical then their hybridized 
modes will repel each other and the crossing will be 
avoided. In our system we expect that the fabricated 
rings would never be perfectly identical, and their differ- 
ent orientation with respect to the waveguide would give 
them slightly different resonant frequencies, thus prevent- 
ing degeneracy. However, our system also has significant 
losses, primarily due to radiation into the forward and 
backward propagating waveguide modes, but also some 
absorption due to metallic and dielectric losses. 

To analyze the more realistic model which includes de- 
tuning and losses, it is convenient to use the theory of 



Morse critical points, which has previously been applied 
to resonators and waveguide^^^'^^oi jj^ j^gf [21] an ap- 
proach was demonstrated which showed the conditions 
for crossing or anti-crossing of modes, and in particu- 
lar demonstrated that losses can counteract an avoided 
crossing. We apply this approach to our system and 
present the key results. For a complete derivation and 
theoretical background, the reader is referred to Ref. 
First we modify the system of equations given by Eqs. ^ 
by introducing a dissipation coefficient F and a detuning 
parameter Sui: 



Qi + 2TQi + (cjo + SluY Qi = -aQ^ 
Q2 + SFQs + K - Q2 = -c^Qi 



P^lQi- (5) 



By taking the time dependence of Qi^2 as exp(jcji), we 
arrive at the following dispersion equation: 



D{uj,e)= (wo + (5w) + j2rw - 



(wo — 5oj) -t- j2ra; — 



(6) 



The solutions of this system are not strictly symmetric 
and antisymmetric, although for small 5uj and for angles 
away from the avoided-crossing they are only slightly per- 
turbed from the original modes. This model neglects the 
differences in radiation losses of the symmetric and anti- 
symmetric modes, which also vary with twist angle. For 
this example we take a detuning parameter 5uj = IO^^wq, 
and a and /3 are given by Eq. Q, with the same coeffi- 
cients used to derive Figs. |4][e,f 

The dispersion curves of the system correspond to 
the condition I?(a;, 9) = 0. In applying the theory of 
Morse critical points, we study the behavior of the func- 
tion D{uj,9) in the neighborhood of the crossing/anti- 
crossing, and do not limit ourselves to values of (w, 9) 
which satisfy the dispersion equation. 

The first step is to find the Morse critical point 
{ujrn,9m) which Satisfies D'^{ujm,9m) = D'g{uJm,9m) = 0. 
In the case of a crossing of modes, this will be the 
point where the dispersion curves cross, in the case of an 
avoided crossing this will be a saddle-point of D{uj,9). 
In general this point must be found numerically, however 
for our system of equations we find a sufficiently accurate 
solution by using perturbation theory to first order in F: 

t^m = ^raO + jT (7) 
WmO = (wq + SUJ^) ^ (8) 

2ao/3iw„o^^o 



cos 



(9) 



In the vicinity of the Morse point we perform a second- 
order Taylor expansion of -D(w, 9) 

D{uj, 9) = e,„) + D'^^{uj - Lo„,f/2 + 

D'^giiJ ~ UJm)i9 - 9rn) + D'^g{9 ~ 9^f /2. (10) 
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For our system the relevant partial derivatives can be 
calculated analytically, but the expressions are long and 
not particularly illuminating. From this local form of the 
dispersion equation, the solutions are given by: 



n" ^_ n" n" 



0mr-2D'M^rn,Orr,) (11) 



This function has a pair of branch points, where the argu- 
ment of the square root becomes zero. Upon substitution 
of the parameters of our system, we find that the branch 
points occur at the following valueJ22l of 9: 



6*61,2 = dm± 2j6ujUJo 



_ .. . (12) 

As shown in Ref. 21, the position of branch points 9bi 2 
on the complex plane determines the behavior of the res- 
onant frequencies of the modes. In brief, if the imaginary 
parts of the two branch points have different signs, then 
the modes undergo an anti-crossing. However, the modes 
cross if the imaginary parts of the 9bi,2 have the same 
sign. This allows us to analyze the eff'ect of detuning 
6uj and losses F on the crossing-anti-crossing behavior 
of the system. To do this, we fix the value of detuning 
5uj = 10~'^a;o and change the losses F. 

We start with weak losses, when F = 1 x IQ^^ujq. This 
case is very close to lossless regime, and the real parts 
of frequencies are plotted in Fig. [sja) , where we see an 
avoided crossing, and in Fig. [sj^b) we see the imaginary 
parts, which are exchanged. 

On the other hand, if we increase losses sufficiently, 
we observe a fundamental change in the behavior of the 
dispersion curves, and the crossing of the real parts of 
the eigenmodes appears. In our system we set F = 
7 X lO'^wo, and the resulting dispersion and absorption 
curves can be seen in Fig. [sf^c) and (d) respectively. If 
we increase F further to 1 x lO^^woi the curves become 
smoother in the vicinity of the Morse critical point, show- 
ing similar character to the numerical data [see Fig. [5je) 
and (f)]. Since to first order uji = F, this case is compa- 
rable to the experimental and numerical results in Fig. [2] 
We have shown how degeneracy (mode crossing) can 
arise in the experimental data, even when the rings are 
not absolutely identical. According to our model, the 
threshold of losses which will mitigate the avoided cross- 
ing (F/cjQ ~ 7 X 10^"^) is similar to the experimentally 
measured values. However this is based on the assump- 
tion that the resonant frequencies of the two rings varied 
by no more than 0.1% {5i^ — lO^'^wo)- The fabrication 
and alignment tolerances in our system are greater than 
this, thus we infer that there is a small avoided crossing 
in our system. It is not possible to measure this directly 
from experiment, since each angle corresponds to a dif- 
ferent sample with different fabrication errors. 

From the numerical results in Fig. and Fig. [5][d), 
we see that even without any fabrication error, the decay 
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FIG. 5. (a) Resonant frequencies and (b) absorption coeffi- 
cients when r = 1 X 10~^aJo, (c) Resonant frequencies and 
(d) absorption coefficients when F = 7 x lO^^oJo, showing 
that losses restore the crossing, (e) Resonant frequencies 
and (f) absorption coefficients when losses are increased to 
r = 2 X 10~^a;o. The projections of the Morse point are 
shown by crosses. 



constants of the hybridized modes are quite different, due 
to different radiation losses. We note that if complete 
degeneracy of both the frequency and decay constants 
of the eigenmodes could be achieved, then the coupling 
between the rings would be completely suppressed. This 
would have important implications for mitigating spatial- 
dispersion, which often inhibits the ability of bulk meta- 
materials to achieve specified properties. 



IV. CONCLUSION 

We have shown in a series of experiments how changing 
the relative rotation between two rings modifies the elec- 
tric and magnetic interaction between them, thus tuning 
the hybridized resonances. Using both numerical sim- 
ulations and an analytical model, that takes into ac- 
count both electric and magnetic interactions between 
the rings, we have shown that there is a crossing where 
the two resonances coexist, at an angle where the elec- 
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trie and magnetic coupling are equal. We have shown 
how experimental errors can cause the crossing of the 
dispersion curve to be avoided, however using the the- 
ory of Morse critical points, we have demonstrated that 
losses can cause the crossing to be restored. 
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